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Turbulent  flows  Involving  chemical  reactions  are  a basic  feature  of  many 
combustion  and  propulsion  systems.  The  development  of  calculation  procedures 
for  turbulent  reacting  flows  requires  the  understanding  and  modeling  of  the 
coupling  between  turbulence  and  combustion.  Second-order  closure  modeling  of 
turbulent  flows  provides  a convenient  framework  for  studying  these  interactions 
between  turbulence  and  chemical  reactions. 

Models  for  the  scalar  probability  density  function  (pdf)  have  to  be  developed 
to  achieve  closure  of  turbulent  transport  equations  for  mixing  and  reacting  flows. 
A delta  function  "typical  eddy"  model  has  been  developed  for  the  joint  pdf  of  the 
scalar  variables.  It  has  been  demonstrated  that  delta  functions  are  a necessary 
part  of  pdf's  in  order  to  attain  the  extremums  of  the  statistical  constraints  on 
the  moments.  The  statistical  bounds  on  a number  of  moments  of  interest  have  been 
derived.  It  has  been  proven  that  a rational  pdf  composed  of  a set  of  delta  func- 
tions alone  can  always  be  constructed  at  any  point  within  the  statistically  valid 
moment  space.  The  model  provides  a good  representation  of  actual  pdf's  in  two- 
species,  variable-density  mixing  flows.  The  model  has  been  directly  compared  to 
experimental  pdf  measurements  and  good  agreement  for  higher-order  moments  has  been 
demonstrated.  It  can  be  shown  that  the  delta  function  pdf  model  is  significantly 
simpler  than  other  proposed  pdf  models  and  Is  more  than  adequate  for  the  closure 
of  the  transport  equations. 
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I.  INTRODUCTION 


Turbulent  flows  involving  mixing  and  reacting  flows  are  a basic  and 
important  feature  of  many  combustion  and  propulsion  devices  such  as  gas  tur- 
bine combustors,  ram  jets,  rocket  engines  and  furnaces.  Any  analytical  model 
of  these  devices  has  to  be  able  to  handle  the  effects  of  the  turbulence  of  the 
flow  and  its  interactions  with  other  physical  and  chemical  processes. 

The  use  of  second-order  closure  modeling  of  the  turbulent  reacting  flow 
holds  great  promise.  Considerable  progress  has  been  made  in  recent  years  in 
the  use  of  second-order  closure  procedures  for  nonreacting  turbulent  flows  and 
improved  predictive  capability  of  these  methods  has  been  demonstrated.  The 
extension  of  this  approach  to  variable-density  mixing  and  reacting  flows  is 
now  feasible,  and  a number  of  modeling  efforts  are  in  progress. 

The  presence  of  finite-rate  chemical  reactions  in  a turbulent  flow  intro- 
duces the  problem  of  proper  modeling  of  many  higher-order  correlations  involv- 
ing scalar  variables  such  as  concentration,  density  and  temperature.  The 
transport  equations  for  the  mean  variables  and  the  second-order  correlations 
are  solved  in  a second-order  closure  procedure  and  these  equations  and  especially 
the  chemical  reaction  source  terms  contain  many  third-order  and  higher-order 
correlations,  such  as  a*8,  p*a,  SB7,  kpa,  etc.  These  correlations  have  to  be 
modeled  in  terms  of  the  lower-order  moments  to  close  the  system  of  transport 
equations.  A convenient  procedure  for  modeling  these  scalar  correlations  is 
to  model  or  calculate  the  probability  density  function  (pdf)  for  the  scalars. 


This  procedure  Is  being  used  at  A.R.A.P.,  as  well  as  by  Rhodes,  et  al.  (1974), 
Bray  and  Moss  (1974),  Lockwood  and  Naguib  (1975),  Libby  (1976),  Bonnlot  (1977a) 
and  Kewley  (1977). 

A number  of  the  above  approaches  require  a number  of  simplifying  assumptions 
and  deal  only  with  one-dimensional  pdf's.  Most  of  them  are  also  restricted  to 
fast  chemistry  and  cannot  handle  finite-rate  chemical  reactions.  The  A.R.A.P. 
model,  called  the  delta  function  "typical  eddy"  model  Is  designed  to  be  appli- 
cable to  finite  rate  chemical  reactions  and  models  the  joint  pdf  for  all  the 
scalars  In  a turbulent  reacting  flow.  The  basic  model  has  been  discussed  by 
Donaldson  (1975)  and  Donaldson  and  Varma  (1976).  Kewley  (1977)  Is  using  the 
A.R.A.P.  model  with  some  modifications  and  has  reported  good  results  on  model- 
ing of  reacting  flows. 

Efforts  are  also  being  made  to  avoid  the  assumptions  regarding  the  shape 
of  the  pdf,  by  directly  solving  transport  equations  for  the  pdf  (Dopazo  and 
O'Brien,  1976,  Bonniot,  1977b  and  Pope,  1976).  However,  these  equations  have 
only  been  solved  for  very  simple  flowfields  and  there  appear  to  be  serious 
modeling  problems  for  some  of  the  terms  In  the  equations.  Spalding  (1977)  has 
recently  proposed  an  alternate  approach  that  Involves  the  calculation  of  the 
pdf  by  following  the  age  history  of  various  eddies. 

This  report  discusses  the  statistical  behavior  of  scalar  variables  in 
turbulent  flows.  It  has  been  shown  that  the  specification  of  lower-order 
moments  leads  to  rigid  constraints  on  the  higher-order  moments.  An  Improved 
version  of  the  "typical  eddy"  model  has  been  developed  that  Incorporates  a 
consistent  procedure  for  the  calculation  of  free  parameters  In  the  model. 

This  eliminates  a major  problem  experienced  with  the  originally  proposed  model. 


The  model  has  been  carefully  tested  against  pdf  measurements  In  a two-species, 
variable-density  mixing  layer  (Konrad,  1976)  and  good  results  have  been  obtained. 
It  Is  claimed  that  a rational,  physically  realistic,  delta  function  pdf  model 
can  be  constructed  at  every  point  In  the  statistically  valid  region  of  the 
moment  space  and  that  this  model  can  be  used  to  predict  third-order  (and  higher- 
order)  moments  to  better  accuracy  than  they  can  be  presently  measured--and  with 
sufficient  accuracy  for  closure  of  the  transport  equations. 

The  model  now  has  to  be  extended  and  tested  for  three  species  reacting  flows 
in  a similar  manner  and  we  anticipate  that  the  model  will  also  be  satisfactory 
for  these  flows. 

The  report  is  orgnized  In  three  sections.  Section  II  derives  the  statis- 
tical bounds  on  the  correlations.  Section  III  details  the  "typical  eddy"  model 
and  the  procedure  for  Its  construction.  The  direct  tests  on  the  model  by  com- 
parison to  experimental  pdf  and  species  intermittency  measurements  are  described 
in  Section  IV. 

This  report  only  briefly  describes  the  mathematical  details  of  the  statis- 
tical analysis  which  is  covered  In  greater  detail  In  a companion  report  (Sandrl, 
et  al.  1978).  For  the  sake  of  completeness,  there  Is  some  unavoidable  duplica- 
tion in  these  two  reports. 


II.  STATISTICAL  CONSTRAINTS  OK  CORRELATIONS 

Basic  statistical  principles  can  be  used  to  obtain  a set  of  constraint  con- 
ditions on  correlations  of  various  variables  in  a turbulent  flowfield.  The 
procedures  for  deriving  these  statistical  constraints  is  discussed  in  this 
section  of  the  report.  The  procedure  Is  quite  general  but  at  the  present  time 
our  main  interest  is  in  scalar  variables  and  more  specifically  in  species  mass 
and/or  molar  concentration  variables. 

The  constraint  conditions  are  useful  in  a number  of  ways.  The  statistical 
constraints  on  first  and  second-order  correlations,  for  example,  a,  ft,  aT7r, 
p'a'  are  important  in  the  question  of  "realizability"  of  second-order  closure 
turbulence  models.  The  lower-order  moments  are  calculated  by  the  solution  of  a 
set  of  modeled  partial  differential  equations.  It  is  possible  that  the  modeled 
equations  may  not  have  solutions  that  are  consistent  with  the  independently 
derived  statistical  constraints,  and  this  will  require  appropriate  corrections 
to  the  models  used  in  the  equations. 

The  statistical  constraints  on  the  third-order  and  higher-order  moments 
are  useful  in  formulating  models  for  these  correlations,  as  will  be  discussed 
in  greater  detail  later  on  in  this  report.  They  are  also  useful  in  determining 
model  sensitivity  and  error-bounds  of  the  modeling  procedure. 

The  statistical  constraints  on  various  moments  of  the  fluctuations  are 
basically  derived  from  conservation  conditions  and  the  Cauchy-Schwarz  Inequality. 
The  Cauchy-Schwarz  inequality  is  a special  case  of  a more  general  inequality. 
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called  Holders'  inequality,  and  in  some  cases  this  more  general  statement  is 
useful . 


Holders  Inequality 

If  f and  g are  arbitrary  positive  variables  and  P is  positive,  then 

Ufp(a)P(a)dcil  f~  j g^ajPtcOdaj 

>.  I f(a)*g(a)P(a 

n 


)da 


provided. 


1.1.1 

p q 


(i) 


(2) 


This  can  be  rewritten  as 


fp)p  M\* 


(3) 


Tne  equality  holds  if  and  only  if. 


fP  - Xg* 


(4) 


Cauchy-Schwarz  Inequality 

This  corresponds  to  p * q * 2 in  the  above  equations. 


2 

T*  g2  > Tg 

The  equality  holds  if  and  only  if, 

f ■ Ag 

Two  other  inequalities  are  useful  in  the  analyses. 


(5) 


(6) 
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Jensen's  Inequality 

If  f(a)  Is  convex  (f“  > 0 throughout  the  domain  of  interest,  0 < a < 1) 


then, 

7T5T>f(S)  (7) 

Tchebytcheff 's  Inequality 

If  f and  g are  similar  functions,  that  is  they  both  increase  or  both  de- 
crease in  the  interval,  0 < a < 1,  then 

7g>?g  (8) 

In  principle,  with  suitable  choices  of  f and  g one  can  derive  all  the 
required  statistical  constraints.  In  practice,  this  is  not  so  easy.  A very 
useful  shortcut  is  the  use  of  a "Renormalization  Theorem".  If  we  have  proven 
that. 


r i 

f 1 1 

f r \ 

f2(A)P1(A)dA 

g2(A)P1 (A)dA  > ( 

fgPj(A)dA 

•'o 

'"'o  ' 

then  it  can  be  proven  that  for  H(A)  > 0, 


r i 

f1  / 

f1  \ 

H f 2 ( A) P] (A)dA 

Hg2(A)P1 (A)dA  > 

HfgP, (A)dA  J 

'o 

^0  ' 

‘■t  / 

(9) 


(10) 


The  proof  is  simple.  Since  Eqn.  9 is  valid  for  all  P^ , one  can  select 
P.j  of  the  special  form  HP^  for  any  P^.  Then  Eqn.  10  follows.  Note  that  the 
reverse  passage  from  Eqn.  10  to  Eqn.  9 is  not  valid  in  general. 

With  the  aid  of  these  mathematical  tools  the  various  statistical  bounds 
for  scalar  variables  of  interest  to  us  can  be  derived.  The  procedures  have  been 
developed  for  a variable-density  two  species  mixing  flow.  Consider  a two  species 
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flow  with  constant  pressure  and  temperature.  Let  a and  6 represent  the  mass 
fractions  and  A and  B,  the  mole  fractions  of  the  two  species.  Assume  a is 
the  heavier  species,  W > W0.  We  define  a normalized  density. 


Vp/RT 


p'  « T'  * 0 


1-f 

a 


0 < A < 1 


The  following  relationships  between  the  variables  can  be  easily  derived 


3 pm  3 pa 


A + 1 


a 3 PaTT-TD 

A complete  second-order  closure  turbulence  analysis  of  a two-species  mixing 
flow  solves  transport  equations  for  the  following  scalar  mean  variables  and 
second-order  moments. 


a,  a2,  p*a,  p*2 


It  can  be  shown  that  these  are  the  only  independent  moments  up  to  second- 
order.  Together  with  the  normal Ization  T * 1,  they  provide  5 independent 
items  of  information  for  construction  of  models  for  higher-order  correlations. 
We  are  interested  in  deriving  the  statistical  bounds  on  the  above  moments  and 
a number  of  the  third-order  moments  that  appear  In  the  transport  equations. 

A simple  conservation  constraint  on  a is 


7 


a + 3 * 1 


or  0 < a 1 

then 

and  a2  < a 

Consider  f = a,  g = 1 in  the  Cauchy-Schwarz  inequality 

_ 2 

a2  1 > a 

Therefore,  the  bounds  on  the  second  moment  a2  are 


Using  the  renormalization  theorem  with  H(cz)  = a,  Eqn.  17  becomes. 


a 

Using  the  renormalization  theorem  with  H(a)  = 1-a,  Eqn.  17  becomes. 


2 


1-a 

Therefore,  the  bounds  on  the  third  moment  a3  for  given  a and  a2  are 


Tl 


The  bounds  on  the  third  moment  given  only  a can  be  derived  using  Jensen's 
Inequality  and  Eqn.  16. 


(22) 


It  Is  an  interesting  exercise  to  compare  the  bounds  on  a5  given  by  the 
two  Eqns.  21  and  22.  The  results  are  shown  In  Figure  1.  The  dotted  lines  are 
the  bounds  on  ot}  for  given  a.  For  given  a,  we  now  |3lck  a2  to  correspond  to 
the  middle  of  Its  allowed  range,  that  is  a2  ■ - g-a-  . Then,  the  solid  lines 
indicate  the  bounds  on  a3  for  specified  a and  a2.  The  important  point  to  be 
noted  is  that  the  bounds  on  a3  when  two  lower-order  moments  are  specified  are 
significantly  narrower  than  the  bounds  when  only  one  lower-order  moment  is 
specified.  This  is  quite  significant  and  leads  to  a very  important  conclusion 
for  our  approach  to  the  modeling  of  the  scalar  probability  density  function. 

This  is  discussed  later  In  detail  in  Section  IV  of  the  report. 

The  statistical  constraints  on  p^a  for  given  a and  a2  are  most  easily 
derived  by  starting  from  the  condition, 

(1 -Aa)2  <_  (1-Aa)  • 1 (23) 

The  lower-bound  is  obtained  by  multiplying  successively  by  p*  and  by  a, 
and  using  the  renormalization  theorem. 


P*(l-Aa) 


a < a-Aa2  p*a 
2 

^ a 

P*a  1 _ --= 
a-Aa2 


(24) 


(25) 
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given  a 
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given  a,  a 


__  _t 
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figure  1.  Comparison  of  the  statistical  bounds  on  a5  for 
specified  lower-order  moments. 
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Multiplying  Eqn.  24  by  B ■ 1-a,  the  upper  bound  on  p*a  can  be  obtained. 


(1-a)  £ (l-a)(l-Aa)  p*(l-a) 

p^T  < J-r  ,A(a-^) 

_TTS  (1-?)  - 


(26) 


Therefore, 


a-Aa^  (1-a)  - A(a-aT) 


(27) 


The  upper  bound  on  p*2  Is  now  considered.  A number  of  useful  relations 
can  be  derived  from  the  equation  of  state, 

P*  - Ap*a  ■ 1 


Then , 


P*aB  » l P*B  - \ B 


P**a  ' 1 P 2 * » P. 


1 

A 


p*2«p  * r p*2p  - I p*e 


P*  2 a2 


, (p*-n 


j 


(28) 


Now, 


T • p*2  > p* 


(29) 


multiplying  by  a and  6 


aB  p*2aB  > p*aB 


(30) 


Substituting  relations  In  Eqn.  28  Into  Eqn.  30,  the  upper  bound  on  p*2  given 
a,  a2  and  p*a  can  be  derived.  The  expression  Is  quite  complicated,  and  Is  not 
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detailed  here.  The  derivation  of  the  lower  bound  of  p*2  for  these  three  given 
moments  does  not  seem  feasible  using  these  manipulations  and  it  has  to  be 
derived  separately  following  the  procedure  discussed  by  Sandri,  et  al.  (1978). 
This  procedure  is  based  on  the  close  relationship  between  the  equality  con- 
dition in  the  Cauchy-Schwarz  theorem  and  the  characteristic  equation  for  a 
discrete  probability  distribution  function. 

The  extremum  values  of  the  statistical  bounds  can  only  be  realized  by 
discrete  probability  density  functions,  that  is,  pdf's  composed  of  delta  func- 
tions. A continuous  function  pdf  must  be  capable  of  going  to  the  limiting 
form  of  delta  functions  in  order  to  be  valid  over  the  entire  statistically 
allowed  range  for  the  moments.  Consider  the  bounds  on  a2.  At  the  upper  bound, 
a2  = a and  the  pdf  can  be  shown  to  correspond  to  two  delta  functions  located 
at  a = 0 and  1.  This  can  be  proven  as  follows: 


a2  = a 


corresponds  to 


f1  f 

•1 

f a2P(a)da  = | 

aP(a)da 

|(a-a2)  P(a)da 

= 0 

*'0 

2 = 0 or  a = 0 and 

a = 1 are  the  only  solutions. 

This 

t ! 

a)  = (1 -a)  6(a) 

+ a 6(a-l ) 

(31) 

as  the  pdf  at  the  upper  bound. 


' V 


P(a)  | 

o 


upper  bound 


At  the  lower  bound. 


lower  bound 


or 


a 


2 


1 

(a2-a  ) P(a)da 


0 


a = a 


and  P(a)  = 6(a-a)  (32) 
is  the  pdf  at  the  lower  bound. 

The  two  limiting  pdf's  are  illustrated  in  the  sketch.  A popular  model  for 
the  pdf  P(a)  is  a Gaussian  or  a clipped  Gaussian  (Lockwood  and  Naguib,  1975). 

For  small  fluctuations,  this  model  does  approach  the  limiting  pdf  correspond- 
ing to  the  lower  bound  on  a2.  However,  it  cannot  attain  the  upper  bound  limit 
corresponding  to  maximal  fluctuations. 

The  same  result  can  be  proven  for  the  upper  and  lower  bounds  of  all  the 
other  moments  of  Interest.  In  fact,  the  bounds  on  the  moments  can  be  easily 
derived  by  considering  the  appropriate  delta  function  pdf.  The  lower  bound 
on  p*2  Is  derived  by  following  this  approach  (Sandrl,  et  al.,  1978).  These 
results  Indicate  that  delta  functions  are  a very  necessary  component  of  the 
pdf's.  Continuous  pdf's  alone  cannot  be  valid  over  the  entire  moment  space. 
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Ideally,  a pdf  composed  of  delta  functions  and  continuous  functions  will  be 
desirable.  However,  our  studies  have  shown  that  rational  pdf's  composed  of 
delta  functions  alone  can  be  constructed  at  every  point  In  the  statistically 
valid  moment  space  and  that  these  pdf's  are  simpler  and  provide  sufficient 
accuracy  for  closure  of  the  transport  equations.  It  Is,  therefore,  unnecessary 
to  develop  the  more  complex  pdf's  and  we  recoronend  the  delta  function  pdf's  be 
used  for  second-order  closure  analyses. 


III.  "TYPICAL  EDDY"  PDF  MODEL 


The  "typical  eddy"  model  is  an  attempt  to  define  the  joint  probability 
density  function  (pdf)  for  all  the  scalar  variables  by  using  the  available 
information  from  a second-order  closure  analysis.  The  pdf  is  represented  by 
a set  of  delta  functions  of  variable  strengths  and  positions  in  the  scalar 
phase  space.  The  strengths  and  positions  of  the  delta  functions  at  every  point 
in  the  flow  are  determined  by  matching  the  model  moments  to  the  values  of  the 
means  and  second-order  correlations  obtained  at  that  point  from  the  solution 
of  their  transport  equations.  The  degrees  of  freedom,  that  is,  the  number  and 
positions  of  the  delta  functions  are  limited  by  the  number  of  independent 
parameters  available  in  the  second-order  closure  analysis. 


Model  Structure  for  Variable-Density  Flows 


In  a two-species,  variable-density  mixing  flow,  the  second-order  closure 
procedure  provides  information  on  four  (4)  independent  moments--a,  a7,  ap*  and 
p*J.  Together  with  the  normalization  condition  that  the  total  probability  must 
be  equal  to  1,  this  provides  us  with  5 pieces  of  information  with  which  to  con- 
struct the  pdf.  After  a large  number  of  tests  we  have  selected  the  following 
model  structure  for  two  species. 
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The  model  has  six  (6)  parameters--e^  thru  e^,  and  <--that  have  to  be 
determined  and  only  5 available  conditions.  Therefore,  the  model  structure 
is  obtained  as  a function  of  the  parameter  o^.  The  model  corresponds  to  the 
following  probability  distribution  function.  The  probability  of  finding  pure 
species  a at  a point  is  given  by  and  the  probability  of  finding  pure  species 
6 is  e Alternatively,  the  strengths  of  the  delta  functions  at  a = 1 and 
a * 0 are  and  e 2 respectively.  The  model  is  composed  of  two  other  delta 
functions  whose  locations  are  variable.  These  delta  functions  of  strengths 
e3  and  e4  are  located  at  a = and  a = tc,  respectively.  The  strengths  of  the 
four  delta  functions  and  the  locations  of  two  of  them  vary  at  different  points 
in  the  flowfield  as  the  values  of  the  means  and  the  second-order  correlations 
are  different.  The  calculation  of  the  model  structure  involves  the  solution 
of  a set  of  nonlinear  algebraic  equations  and  in  the  past  has  involved  a great 
deal  of  effort.  We  have  now  developed  a very  convenient  procedure  for  the  solu- 
tion of  the  model  structure  that  exploits  the  mathematical  properties  of  the 
delta  functions.  The  pdf  for  the  model  structure  can  be  written  as, 

P(a)  * S(a-l)  + e ^ 6(a)  + 6(a-a3)  + 6(a-tc)  (33) 

ftl 
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It  is  convenient  to  introduce  two  delta  functions  in  density  space  in  the 
above  expression.  Using  the  equation  of  state  (Eqn.  13),  the  pdf  can  also  be 
written  as, 

P(a)  - c1  6(p*-  + e2  6(p*-1)  + e3  <5(a-a3)  + e4  6( a-tc)  (34) 

With  repeated  use  of  Dirac's  identity,  that  is, 

(a-a)  6(a-a)  = 0 


the  following  equation,  called  the  characteristic  equation  of  the  pdf  is  ob- 
tained : 


(p*-  77^)  (p*-1  ) (u-a3)  (a-tc)  * 0 


(35) 


Taking  the  time  average  of  this  relation  gives  an  expression  for  k. 

<or  (p*-P*(  1+  )>  - a3  <a(p£-p*(l+  ) + yr£  )> 

<7  (:;-?*(!+  777  ) + T7£  )>  - a3  < ( P*P*  ( 1 + TTa  } + TTa 


(36) 


With  the  use  of  a number  of  relations  derived  from  the  equation  of  state 
(see  Eqn.  28),  the  above  expression  can  be  simplified  to  express  < as  a func- 
tion only  of  the  four  moments  a,  a2,  p*a,  p*  and  the  parameter  a^. 

The  cell  sizes  thru  e4  can  be  calculated  as  follows.  Multiply  Eqn.  34 
by  (p*-l),  (0-03)  and  (a-tc).  Then, 

(p*-l  ) (a-a3)  (a-tc)P(a)  * e ^ ( p„- 1 ) (a-a3 ) (a-< ) <5  ( p*-  ^7^  ) 

* E1  (77^  ' 1 )0“«3)0-k)S(p*-  77^) 
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Integrating  over  a from  0 to  1, 


<(p*-l  ) (a-a3)  (a-<)>  = ei  " 1 ) ( 1 -<*3)  ( 1 -k) 

<(p*-l  )(a-a3)(a-ic)> 

^ (jTh  “1 ) ( 1 -a3)  ( 1 “K ) 


The  expressions  for  the  other  e's  can  be  derived  in  a similar  manner. 

<(p*-  y7^)(a-a3) (q-k)> 

e2  = 7\  JZ\ 

a3  K 1 - A 

<(p*-1)(p*-  ( a-K  ) •” 

3 ^3'k)(u1^  • I^A) 

<(p*-1)(p*-  (a_a3)> 

4 (K_a3)(,r^ 


With  a little  bit  of  algebra  and  the  use  of  Eqn.  28,  the  above  expression 
for  the  cell  sizes  can  also  be  written  as  a function  of  the  four  independent 
moments  and  the  parameter  a-j. 

A computer  program  has  been  written  to  calculate  the  complete  two  species 
model  structure  for  given  values  of  the  moments.  The  program  scans  the  entire 
statistically  valid  moment  space,  and  calculates  the  cell  sizes  and  < for 
values  of  e3  from  0 to  1.  A solution  for  the  model  pdf  is  physically  correct 
only  if  all  the  cell  sizes  (strengths  of  the  delta  functions)  are  positive  and 
the  location  ic  is  between  0 and  1.  Valid  solutions  of  the  model  structure  have 
been  obtained  at  all  points  within  the  statistically  valid  region.  In  general. 
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these  valid  solutions  are  found  for  a limited  range  of  values  of  the  parameter 
aj.  An  arbitrary  specification  of  (for  example,  ^ or  a)  is  not  valid  over 
the  entire  range  of  moments.  This  is  consistent  with  our  experience  with  an 
earlier  version  of  the  "typical  eddy"  model  for  two  and  three  species  flows 
(Donaldson  and  Varma,  1976)  in  which  the  proportions  in  the  mixed  cells  were 
specified  to  be  At  some  points  in  the  flow,  this  led  to  the  occurrence  of 
negative  probabilities  and  unphysical  models. 

A reasonable,  consistent  criterion  is  still  necessary  for  the  selection  of 
dj.  Of  course,  if  one  additional  piece  of  information  was  available,  would 
be  uniquely  determined.  However,  in  an  analysis  limited  to  the  second-order 
closure  level,  no  other  independent  information  is  available  and  a reasonable 
assumption  has  to  be  made,  can  be  picked  at  any  point  within  its  valid 
range  of  values.  The  results  of  the  computer  program  calculations  prove  that 
at  any  point  in  moment  space,  there  are  very  small  changes  in  quantities  of 
interest,  for  example,  third-order  moments,  over  the  entire  valid  range  of 
values  of  a^.  A few  examples  of  these  results  are  discussed  later  in  this 
section.  It  is,  therefore,  concluded  that  the  exact  choice  of  is  not  very 
critical . After  careful  study  and  comparison  of  the  results  of  various  assump- 
tions with  experimental  measurements  of  the  pdf  (Section  IV)  we  have  adopted 
the  following  empirical  assumption  for  selection  of  o^. 

"The  parameter  is  selected  within  its  valid  range 
at  the  point  of  minimum  mixedness  or  minimum  entropy  of 
mixing." 

The  transport  equations  for  the  means  and  the  second-order  correlations 
incorporate  the  entropy-increasing  effects  of  molecular  transport  processes  and 


I 


9 


19 


the  values  of  the  correlations  at  a point  Include  this  entropy  change.  The 
entropy  of  mixing  for  the  pdf  is  calculated  as, 


n - ^ Ci  [f  Pi^l  iT  10  P1ai  ' P1*1  ^ 10  Pi6lJ  (41) 

The  entropy  of  mixing  n is  a function  of  the  parameter  o^.  The  true 
entropy  of  the  flow  at  a point  must  range  between  the  minimum  and  maximum 
values  obtained  over  the  valid  range  of  a^.  However,  the  true  entropy  of 
mixing  (which  can  be  considered  to  be  a function  of  a set  of  higher-order 
moments  of  the  pdf)  is  unknown  in  a second-order  closure  analysis  and  we 
make  the  choice  of  to  correspond  to  the  point  of  minimum  entropy  of  mixing. 

This  ensures  that  the  model  entropy  is  never  larger  than  the  true  entropy  of 
the  system.  The  modeling  procedure,  thereby,  does  not  add  entropy  to  the 
system. 

The  concept  of  minimum  entropy  of  mixing  is  the  same  as  minimum  mixedness. 
This  is  easily  illustrated  by  considering  the  following  simple  pdf's.  For 
specified  a the  entropy  of  mixing  (Eqn.  41)  is  always  zero  for  the  pdf  a, 
while  it  is  the  maximum  for  pdf  b.  The  pdf  a correspo*  ?s  to  minimum  mixedness. 


that  is,  there  is  zero  probability  of  finding  molecularly  mixed  regions  of  a 
and  8.  Pdf  b corresponds  to  the  most  mixed  state. 

The  selection  of  at  the  position  of  minimum  mixedness  is  also  supported 
by  the  observed  large-scale  eddy  structures  in  turbulent  flows.  Consider  the 
mixing  of  two  initially  separate  streams  of  a and  8.  The  large  eddies  in  the 
flow  cause  stretching  and  folding  and  intertwining  of  the  streams,  but  do  not 
cause  an  increase  in  the  entropy  of  the  system.  The  entropy  of  mixing  increases 
due  only  to  the  molecular  mixing  processes  at  the  interfaces  of  the  two  streams. 
Selecting  the  minimum  entropy  model  is  consistent  with  stressing  the  dominance 
of  the  large-scale  structures  in  the  flow. 

Ultimately,  however,  the  choice  of  must  be  considered  to  be  an  empirical 
choice.  An  alternative  selection  criteria  that  has  been  tested  to  some  extent 
is  to  use  the  middle  value  in  the  allowed  range  of  a 3 at  every  point.  This 
model  is  able  to  predict  third-order  moments  better  than  the  minimum  entropy 
model,  but  other  comparisons  of  model  predictions  and  experiments  indicate 
that  the  minimum  entropy  model  shows  the  best  overall  agreement  with  the  measure- 
ments for  the  one  flow  we  have  studied  in  detail  (Section  IV).  Further,  the 
midrange  of  model  involves  greater  complexity  in  model  construction  than 
the  minimum  entropy  model  and  the  latter  has  been  selected  as  the  model  to  be 
used  in  our  second-order  closure  program.  It  has  also  been  shown  that  third- 
order  moments  are  not  very  sensitive  to  the  choice  of  in  our  pdf  model.  This 
is  because,  given  4 lower  order  moments,  the  third-order  moments  are  tightly 
constrained  and  as  long  as  the  model  is  statistically  valid,  the  third-order 
moments  can  be  predicted  with  good  accuracy. 
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Examples 

Two  examples  of  the  results  obtained  from  the  computer  program  are  now 
discussed.  Consider  a flow  of  two  gases  a and  B of  molecular  weights  Wq  = 32 
and  WQ  = 2.  Consider  a point  in  the  statistically  valid  moment  space. 


a 

a2 

P a 

pI 

0.2 

0.1 

1.0 

12.0 

The  valid  range  of  is  found  to  be  0.12  to  0.27  from  the  computer  pro- 
gram. Figure  2 shows  the  results  for  the  third  moments  a3,  p*3  and  the  entropy 
of  mixing  as  a function  of  a^.  The  figure  shows  that  the  third-moments  are 
nearly  constant  over  the  valid  range  of  a^. 

Consider  another  point  in  the  moment  space  where  the  four  independent 
moments  have  the  following  values. 


a 

a2 

P a 
r* 

IV 

0.6 

0.4 

2.0 

12.0 

The  valid  range  of  is  now  found  to  be  0.37  to  0.61.  The  results  for 
the  moments  and  the  entropy  of  mixing  are  shown  in  Figure  3,  and  again  demon- 
strate the  insensitivity  of  the  higher-order  moments  to  the  selection  of  a3 
within  the  stylistically  valid  region. 
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Model  Structure  for  Constant-Density  Flows 

In  a constant-density  flow,  the  second-order  closure  procedure  provides 
information  on  only  two  (2)  Independent  moments--a  and  a2--and  the  following 
model  structure  is  used. 


0 €, 


cjj  is  again  a free  parameter.  The  model  structure  solution  involves  only 


linear  algebraic  equations  in  this  case.  Namely, 


C1  * c2  * e3  ’ 1 


1 ' e 1 + °-c2  * °3 ’ e3  " a 


1 •c1  + 0-e2  + * e 3 = a2 


These  can  be  easily  solved  and  lead  to, 


a - aa, 


. , — a - a‘ 

e?  s 1 - a - 

1 1 - a 


. . a - a‘ 


a3  is  selected  at  the  position  of  minimum  entropy  of  mixing. 
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IV.  COMPARISON  OF  "TYPICAL  EDDY"  MODEL  WITH  PDF  MEASUREMENTS 


Normally  a turbulence  model  of  this  type  Is  tested  within  the  framework 
of  a complete  solution  of  the  second-order  closure  transport  equations  and  com- 
parison of  the  results  of  the  calculation  against  experimental  measurements. 

This  Is  necessary  because  detailed  measurements  of  turbulence  correlations  are 
usually  not  available.  However,  it  has  been  possible  to  make  a direct  test  of 
the  "typical  eddy"  pdf  model.  This  has  enabled  us  to  test  the  adequacy  of  the 
model  and  to  ascertain  its  sensitivity  to  various  parameters.  The  comparisons 
have  enabled  us  to  provisionally  select  the  minimum  entropy  of  mixing  criteria 
for  the  selection  of  the  model  parameter  a^.  The  direct  tests  of  the  model 
have  greatly  increased  our  understanding  of  the  model  structure. 

Konrad  (1976)  has  made  a series  of  detailed  measurements  of  the  pdf  in 
three  shear  layer  flowfields.  The  measurements  have  been  made  in  (i)  a He-N^ 
shear  layer  with  a velocity  ratio  ^/u-i  s 0-38  and  a density  ratio  P2|/pl  = ?.0. 
(11)  a uniform  density  He+Ar-N^  shear  layer  with  velocity  ratio  of  u?/u1  = 0.38 
and  (iii)  a wake  flow  in  which  the  gases  have  equal  free-stream  velocities  and 
a density  ratio  of  7.  All  the  flows  are  low  speed  and  two-dimensional.  The 
measurements  consist  of  mean  velocity,  mean  concentration,  concentration 
fluctuation  correlations,  species  intermittency  functions  and  the  species  pdf's. 
The  latter  two  parameters  are  used  for  the  direct  model  comparison. 

The  pdf  measurements  are  used  to  determine  the  values  of  the  four  lower- 
order  moments--a,  a2,  ap*  and  p^--that  are  needed  to  construct  the  two-species 
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variable-density  "typical  eddy"  model.  The  measured  pdf  Is  also  used  to  calcu- 
late experimental  values  of  a number  of  third-order  correlations  of  Interest, 
such  as  aB2,  AB2  and  A2B.  The  third-order  moments  are  also  computed  from  the 
model  and  the  results  compared  to  the  measurements.  The  valid  range  of  the 
free  parameter  a^  corresponds  to  the  upper  and  lower  bounds  on  the  third  moments. 
The  use  of  the  constant-density  "typical  eddy"  model  In  a variable-density  flow 
Is  also  Investigated.  In  this  case  the  model  Is  constructed  by  only  using  a 
and  a2. 

The  data  for  the  He-N2  shear  layer  with  a velocity  ratio  of  0.38  and  a 
density  ratio  of  7 has  been  most  extensively  used  for  model  testing.  Figure  4 
shows  the  results  for  the  third  moment  aB2.  The  ordinate  of  the  figure  Is 
“^model  expt'  The  absc^ssa  the  normalized  coordinate  across  the  shear 
layer.  The  dotted  lines  show  the  upper  and  lower  bounds  on  the  third  moment 
when  only  two  lower-order  moments  are  used  for  the  model  construction.  In  this 
case  0tP2mo(jei  has  a large  range  of  possible  values  and  some  models  within  the 
statistical  range  can  lead  to  significant  errors  compared  to  the  experiments. 
However,  when  4 lower-order  moments  are  specified,  the  model  values  of  aB2  are 
tightly  constrained  as  shown  by  the  solid  lines.  In  fact,  now  any  statistically 
valid  model  Is  able  to  calculate  aB2  (and  other  third-order  moments)  to  better 
than  10*  accuracy.  This  Is  significantly  better  accuracy  than  tt.e  expected 
error-bounds  on  experimental  measurements  of  third-order  moments. 

For  this  particular  flowfleld,  the  lower  bound  (solid  line)  on  aB2  corres- 
ponds to  the  minimum  entropy  of  mixing  anr  the  upper  bound  (solid  line)  corres- 
ponds to  the  maximum  entropy  of  mixing.  An  alternate  empirical  choice  for 
In  the  model  Is  suggested  by  Figure  4;  namely,  to  select  a^  at  the  middle  of 
Its  allowed  range.  The  results  for  this  choice  of  a3  are  also  shown  In  the 
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figure.  This  leads  to  the  most  accurate  results,  with  the  third  moment  pre- 
dicted to  better  than  2%  accuracy  across  the  entire  flowfield.  However,  other 
comparisons  with  data  suggest  better  results  with  the  minimum  entropy  of 
mixing  criteria.  Further,  the  selection  of  the  midrange  value  of  a3  Involves 
additional  complexity  In  the  computer  program  and  for  the  present  we  are  pro- 
posing to  use  the  minimum  entropy  criteria.  Note  that  this  assumption  still 
enables  us  to  predict  aB2  to  within  10%  of  Its  measured  value  for  a variable 
density  flow.  This  Is  quite  adequate.  The  use  of  the  midrange  value  of  a3 
is  an  option  available  to  us  If  further  comparison  with  measurements  In  other 
flows  suggest  the  need  for  a more  accurate  selection  of  a3. 

Figure  5 presents  the  results  for  the  third  moment  A2B  In  the  He-N^  shear 

layer.  The  variable-density  "typical  eddy"  model  with  four  moments  specified 
Is  able  to  predict  third-order  moments  to  better  than  10%  accuracy.  The  lower 
bound  solid  line  corresponds  to  the  minimum  entropy  choice  for  a,. 

The  results  for  the  variable  density  He-N2  wake  flow  (u2/u1  ■ 1.0, 

P2/Pl  = 7.0)  also  lead  to  the  same  conclusions.  Figure  6 shows  the  model  pre- 
dictions for  aB2  compared  to  the  experimental  measurements. 

A qualitative  comparison  of  the  "typical  eddy"  delta  function  pdf's  and 
the  measured  pdf's  Is  shown  In  Figures  7 and  8.  Figure  7 shows  the  results  for 
the  He-N2  shear  flow  with  a3  selected  on  the  basis  of  minimum  entropy  of  mixing. 
The  experimental  pdf  Is  reproduced  from  the  report  by  Konrad  (1976).  The  delta 
functions  In  the  model  are  shown  at  their  correct  locations  in  the  concentration 

space  and  the  relative  strengths  of  the  delta  functions  are  correct.  The  absolute 

magnitude  of  the  strength  of  the  delta  functions  cannot  be  directly  plotted  on 
this  graph  and  we  have  selected  a convenient  scale  for  depicting  the  delta 
functions.  The  delta  function  representation  of  the  pdf  seems  to  be  capable 
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Figure  5.  Comparison  of  model  predictions  and  experiments  for  the 
third  moment  TRL  He-N2  shear  layer. 
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of  capturing  the  important  features  of  the  experimental  pdf  structure. 

Figure  8 shows  the  corresponding  results  when  a ^ is  selected  on  the  basis 
of  the  maximum  entropy  of  mixing.  A comparison  of  Figures  7 and  8 suggests 
that  the  minimum  entropy  criteria  leads  to  better  results  than  the  maximum 
entropy  criteria. 

The  results  for  the  third-order  moment  predictions  and  the  pdf  represen- 
tations for  the  uniform-density  He  + Ar  - Ng  shear  layer  ^/u^  = 0.38)  are 
now  discussed.  In  this  case  only  two  lower-order  moments  a,  a2  are  available 
to  construct  the  model  and  the  bounds  on  aP2  are  considerably  wider  as  shown 
in  Figure  9.  The  model  predictions  for  the  minimum  entropy  of  mixing,  maximum 
entropy  of  mixing  and  setting  = a are  shown  in  the  figure.  The  last  model 
is  of  historical  interest  to  us.  The  accuracy  in  the  prediction  of  a£2  is  not 

very  good  with  any  of  the  models,  but  the  minimum  entropy  model  has  the  least 

error.  Other  comparisons  with  data  discussed  later,  also  support  the  minimum 
entropy  model  for  this  flowfield.  The  midrange  value  of  ci^  has  not  yet  been 
investigated  in  this  flow  but  is  expected  to  yield  the  best  results  as  far  as 
third-moment  predictions  are  concerned.  The  delta  function  pdf's  for  three 
choices  of  are  compared  to  experimental  measurements  in  Figures  10,  11,  and 
12.  The  maximum  entropy  model  is  rather  poor  in  simulating  the  large  probabil- 
ities of  pure  species  at  the  edges  of  the  flowfield.  The  minimum  entropy  model 
Is  qualitatively  somewhat  better  than  the  = a model. 

The  final  model  comparison  with  the  experiments  is  for  the  probability  of 
finding  pure  species  at  various  points  across  the  flowfield.  Figure  13  pre- 
sents the  results  for  the  uniform  density  shear  layer  flow.  The  predictions 

for  = a model,  maximum  entropy  and  the  minimum  entropy  of  mixing  model  are 
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Figure  9.  Comparison  of  model  predictions  and  experiments  for  the 
third  moment  aB7.  Constant-density  He+Ar-N2  shear  layer. 


Figure  10.  Comparison  of  delta  function  pdf  with  measured  pdf.  Minimum  entropy  of  mixing 
model.  Constant-density  He+Ar-N2  shear  layer. 
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Flqure  13.  Comparison  of  model  predictions  with  experiments 
9 for  the  probability  of  finding  pure  species. 

Constant-density  He+Ar-N2  shear  layer. 
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compared  to  the  data.  The  theoretical  predictions  use  a cutoff  value  of  0.9, 
that  is,  a cell  In  the  model  is  considered  to  be  a pure  cell  of  species  a if 
the  proportion  of  a in  it  is  greater  than  or  equal  to  0.9.  This  is  an  arbi- 
trary, but  reasonable,  criterion. 

Figure  13  shows  that  the  maximum  entropy  model  is  not  at  all  satisfactory, 
while  the  results  for  the  two  other  models  are  comparable.  The  = a model 
was  used  in  the  past,  and  although  it  is  satsifactory  for  constant  density  flows, 
it  fails  in  the  variable  density  cases  and  is  not  consistently  valid.  The  mini- 
mum entropy  of  mixing  selection  criterion  for  gives  good  results  for  this 
constant-density  flow. 

Figures  14  and  15  present  the  results  for  the  pure  species  probabilities 
in  the  He-N^  variable-density  shear  layer  flow.  The  results  in  Figure  14 
using  the  minimum  entropy  of  mixing  for  selection  of  appear  to  be  signifi- 
cantly better  than  the  use  of  the  maximum  entropy  or  the  midrange  value  of  a^. 

This  limited  comparison  with  detailed  pdf  measurements  indicates  that  the 
empirical  model  of  selecting  to  correspond  to  the  minimum  entropy  of  mixing 
shows  the  best  overall  agreement  with  the  data.  An  alternative  model  of  selec- 
ting aj  at  its  midrange  value  is  better  in  some  aspects,  but  at  this  time  the 
increased  complexity  of  its  calculation  is  not  worthwhile.  Comparison  with 
other  measurements  for  different  flow  conditions  is  desirable  for  further  model 


evaluation,  but  the  current  model  provides  sufficient  accuracy  for  the  predic- 
tion of  third-order  moments  and  the  closure  of  the  transport  equations. 
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Figure  15.  Comparison  of  model  predictions  with  experiments  for 
the  probability  of  finding  pure  species.  Maximum 
entropy  of  mixing  and  a.  at  midrange  models. 

He-N2  shear  layer. 


V.  CONCLUSIONS 


The  direct  comparison  of  the  4 cell  delta  function  "typical  eddy"  model 
with  pdf  measurements  has  demonstrated  that  this  simple  model  is  quite  adequate 
for  achieving  closure  of  the  transport  equations.  For  a two  species,  variable 
density  mixing  flow  we  have  shown  that, 

• Delta  functions  are  a necessary  part  of  the  pdf  in 
order  to  attain  the  extremum  values  of  the  statis- 
tically valid  moment  space.  Continuous  pdf's  alone 
(without  delta  functions)  are  not  sufficient,  and 
can  be  statistically  invalid. 

• A rational  pdf  composed  of  a set  of  delta  functions 
can  always  be  constructed  anywhere  within  the  statis- 
tically valid  moment  space. 

• A pdf  composed  of  delta  functions  alone,  that  is,  the 
A.R.A.P.  "typical  eddy"  model,  when  using  all  the 
available  information  in  a second-order  closure  calcu- 


lation, can  predict  higher-order  correlations  to  an 
accuracy  better  than  they  can  be  measured.  It  is, 
therefore,  unnecessary  to  construct  more  complex, 
continuous  pdf's  for  the  purpose  of  achieving  closure 
of  the  transport  equations. 


The  model  now  has  to  be  extended  to  reacting  flows  Involving  three  and 
more  species.  The  mathematical  framework  that  has  been  established  In  the 
development  of  the  two  species  model  will  be  extensively  used  In  this  effort. 
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developed  to  achieve  closure  of  turbulent  transport  equations  for  mixing  and 
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rational  pdf  composed  of  a set  of  delta  functions  alone  can  always  be  con- 
structed at  any  point  within  the  statistically  valid  moment  space.  The  model 
provides  a good  representation  of  actual  pdf’s  in  two-species,  variable- 
density  mixing  flows.  The  model  has  been  directly  compared  to  experimental 
pdf  measurements  and  good  agreement  for  higher-order  moments  has  been 
demonstrated.  It  can  be  shown  that  the  delta  function  pdf  model  is  signifi- 
cantly simpler  than  other  proposed  pdf  models  and  is  more  than  adequate  for 
the  closure  of  the  transport  equations.^ 
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